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The standard technique for measuring the phase of a single-mode field is heterodyne detection. 
Such a measurement may have an uncertainty far above the intrinsic quantum phase uncertainty 
of the state. Recently it has been shown [H. M. Wiseman and R. B. Killip, Phys. Rev. A 57, 2169 
(1998)] that an adaptive technique introduces far less excess noise. Here we quantify this difference 
by an exact numerical calculation of the minimum measured phase variance for the various schemes, 
optimized over states with a fixed mean photon number. We also analytically derive the asymptotics 
for these variances. For the case of heterodyne detection our results disagree with the power law 
claimed by D'Ariano and Paris [Phys. Rev. A 49, 3022 (1994)]. 
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I. INTRODUCTION 

It is well known that it is not possible to make 
quantum-limited measurements of the phase of an elec- 
tromagnetic field using linear optics and photodetectors 
jij. The standard method of making a (non-quantum- 
limited) phase measurement is by heterodyne detection. 
Heterodyne detection involves combining the field to be 
measured with a much stronger local oscillator field which 
has a frequency detuned by a small amount A, and mea- 
suring the intensity of the resultant field. In a typi- 
cal experimental implementation, the mode to be mea- 
sured is passed through a 50/50 beam splitter, in order 
to combine it with the local oscillator. The difference 
photocurrent from the two output ports of the beam 
splitter yields a measurement of the phase quadrature 
X<j> = ae~ J * + a^e 1 *, where a is the annihilation opera- 
tor of the mode to be measured and <& is the phase of the 
local oscillator. In heterodyne measurement the detuning 
A is typically chosen to be large enough that the phase of 
the local oscillator cycles many times over the course of 
the measurement, so as to measure all quadratures with 
equal accuracy. 

In a heterodyne phase measurement of a coherent state 
of amplitude a, the variance in the measured phase 4> 
scales as [0 
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This is twice the intrinsic uncertainty of the coherent 
state. An improved phase measurement can be made if 
one has an initial estimate of the phase. Then one would 
choose the local oscillator phase $ to be equal to ip+ir/2, 
where ip is the initial estimated phase of the field to be 
measured. This is known as homodyne detection. If the 
phase is unknown before the measurement, then one can 
still apply this idea by adjusting the phase of the local 
oscillator during the course of the measurement based on 
an estimate of the phase from the measurement results 
so far. 



These adaptive phase measurements have been dis- 
cussed in a series of papers by Wiseman || and Wiseman 
and Killip |], [§. In Ref. § it was shown that via adap- 
tive measurements of a coherent state one can obtain a 
phase uncertainty of 



(1.2) 



Here we have taken the coherent amplitude a to be real. 
The first term above is the intrinsic uncertainty of the 
coherent state, and the second term is the extra phase 
uncertainty introduced by the phase measurement. This 
result implies that even though there is not much im- 
provement in using adaptive phase measurements on co- 
herent states, if the input state has reduced phase uncer- 
tainty then adaptive measurements will produce a phase 
measurement with far less uncertainty than a heterodyne 
measurement. 

To better quantify the improvement offered by adap- 
tive measurements, it is necessary to consider the vari- 
ance of states that have been optimized for minimum 
phase variance under various measurement schemes. 
There must be a constraint placed on optimizing states, 
otherwise a phase eigenstate will be obtained. The two 
main ways in which to constrain the states are by trun- 
cating the photon number and by fixing the mean photon 
number. The case of truncated photon number was con- 
sidered analytically in Ref. and numerically in Ref. 
jjj. The case of a fixed mean photon number is deter- 
mined analytically in Sec. IV below. 



General optimized states are difficult to work with, as 
they cannot be generated experimentally. Also, in nu- 
merical integration of the stochastic differential equations 
arising from the various measurement schemes (to be con- 
sidered in future work), calculations must be performed 
on the entire state for general optimized states. This 
becomes prohibitively time consuming for large photon 
numbers. Squeezed states are far more practical, as they 
are routinely generated experimentally, and in numeri- 
cal integration only the two squeezing parameters need 
be considered. The theory for optimized squeezed states 
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is considered in Sec. [v| The analytical results for both 
general states and squeezed states are tested num erica lly 
in Sec. VI, and the implications discussed in Sec. VII. 



II. POMS AND PHASE MEASUREMENTS 

In quantum-mechanical systems, the most general way 
of obtaining the probability of some measurement result 
E is by the expectation value of an operator F(E), i.e., 



P(E) = Tr[pF(E)], 



(2.1) 



where p is the state matrix for the system. If the set 
of all possible measurement results is 57, it is evident 
that P(£l) = 1 for all p, which implies that F(il) = 1. 
Thus F(E) can be called a probability operator, and the 
mapping E i— > F defines a probability operator measure 
(POM) on |. 

For phase measurements on a single-mode field, the 
general form of the POM is J2j 



1 



OO 



Ztt * — ' 

n,rn— 



m)(n\e l 



(2.2) 



where if is a positive-semidefinite symmetric matrix with 
all entries positive, and \m) is a number state of the field. 
In this case fi = [0, 2ir) and the completeness relation is 



\ d<j>F{<f>) =diag[H o, H 11 ,H 22 , . 
Jo 



.] = 1. (2.3) 



For ideal phase measurements all elements of the H ma- 
trix are equal to 1, whereas for physical measurements 
the off-diagonal elements will generally be less than 1 || . 

The accuracy of phase measurements can be quanti- 
fied in a number of ways Q which agree in the limit 
of small phase variance provided the phase distribution 
is narrowly peaked (as it will be in the examples we con- 
sider). When the mean phase is zero it is easiest to define 
the phase variance as 

i<f>\ 



V(</>) =2-2 (cos(/>) = 2(1 - Re(e^)) 



(2.4) 

If we evaluate (e 1 ^) using Eq. ( |2.2| ) for an arbitrary pure 
quantum state \tp), we obtain 



= j^(^|m)<m + l|^)B» 
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Therefore, the phase variance only depends on the off- 
diagonal elements H m ^ m+ i, and we may characterize a 
phase measurement by the vector 



h(m) 
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For most phase measurements, we have for large photon 
numbers H 



h(m) 
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Then for states with a reasonably well-defined mean pho- 
ton number n the total phase variance is given approxi- 
mately by §| 



V((j)) PS ^intrinsic (</>) + 2h(n) 



(2.8) 



where Vi n trinsic(</ ) ) is the phase variance which would re- 
sult from an ideal or canonical measurement of the phase 
H pi- Evidently h(n) is a measure of the excess phase 
noise introduced by the measurement, and it would be 
desirable to make this as small as possible. 



III. ADAPTIVE MEASUREMENTS 

In a real optical experiment one cannot directly mea- 
sure the phase, but would rather estimate it from a pho- 
tocurrent record. Here the photocurrent would be de- 
rived by combining the mode to be measured with a 
local oscillator via a 50/50 beam splitter. Such mea- 
surements have been called dyne measurements, as they 
include homodyne and heterodyne as well as adaptive 
measurements ]lj . The signal of interest is the difference 
between the photocurrents at the two ports, which we 
define by 



lit) = lim lim 

St — >0 (8 — >og 



SN+(t) -SN-(t) 



(3.1) 



where SN± are the increments in the photocounts at the 
two detectors in the interval [t, t + 5) and f3 is the ampli- 
tude of the local oscillator. 

Say the mode to be measured has an (assumed posi- 
tive) envelope u{i) normalized such that J u{t)dt = 1. 
For simplicity, define a scaled time v by 



u(s)ds, 



(3.2) 



so that v € [0, 1). Then it turns out that at time v there 
are two sufficient statistics for the measurement record: 

A v = [ I(u)e^ u) du, B v = — f e 2l * (u) du. (3.3) 
Jo Jo 

These are sufficient statistics in the sense that the POM 
for the measurement is a function of A v and B v only |l0| . 
This means that the best estimate for the phase at time 
v need depend on the measurement record only through 
the complex numbers A v and B v . 

Consider the simple case where the mode is initially in 
a co herent state \a). We can then take the limits in Eq. 
(O) to obtain 



I(v) =2Re(ae- l * w ) +£(«), 



(3.4) 



where £(u) = dW(v)/dv is the shot noise ||. From this 
we can evaluate A v as 

A v = av-a*B v + f e 2l * (tl) >dW(u). (3.5) 
Jo 
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If one ignores the final term (which has zero expectation 
value), it is easy to see that arga = &rg(vA v + B V A*). 
Other arguments Q suggest that this is generally the 
best estimate for the phase at time v. If B v is small (as 
it is in heterodyne detection), this can be approximated 
by sxgA v . The adaptive phase measurements that were 
analyzed in Refs. g] and Q use axgA v as the phase es- 
timate during the measurement, setting 



$(«) = argA„ +7I-/2. 



(3.6) 



The main motivations for this choice are that it gives a 
feedback algorithm which would be easy to implement 
experimentally, and that it is mathematically tractable. 

If the phase estimate aigA is used at the end of the 
adaptive measurement, the resultant phase measurement 
is actually worse than a heterodyne phase measurement. 
This phase measurement is called an adaptive mark I 
phase measurement, and it was found j^j that 

~ o 179 ' ( 37 ) 
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which compares with 
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If the phase estimate axg(A + BA*) is used at the end 
of the phase measurement, a far better phase measure- 
ment is obtained. This phase measurement is called an 
adaptive mark II phase measurement, and yields g] 

hu{m) ~ Tei^' (3 - 9) 

which is considerably better than the standard (hetero- 
dyne) result. However, as noted in Sec. Q, this dramatic 
improvement can only be seen if one starts with a state 
having intrinsic phase variance much less than that of a 
coherent state, for which intrinsic = 2ft.hot(a 2 ). 

In practice, imperfections in the equipment can intro- 
duce phase uncertainties that scale as TT 1 |J. The main 
problems are inefficient detectors and time delays. From 
Rcf . , the extra phase variance introduced by inefficient 
detectors is 



Arjn 



(3.10) 



where 77 is the efficiency. Again from Ref. |2j the extra 
phase variance introduced by a time delay Sv is 



8v 

SVn(<i>) ~ -=-■ 
n 



(3.11) 



This means that when the photon number is sufficiently 
large, above (1 — ?/)~ 2 or (<ki)~ 2 , the introduced phase 
variance will always scale as TT 1 . Nevertheless, adap- 
tive measurements will always give more accurate results 
than heterodyne measurements, and it is of interest to 
know how accurate phase measurements can be made in 
principle, since future technological advances may greatly 
reduce the problems of detector inefficiencies and time 
delays. 



IV. GENERAL OPTIMIZED STATES 

The fairest way to compare different phase measure- 
ment schemes is to consider the phase variances for states 
optimized for minimum total phase variance for the par- 
ticular measurement under consideration. The two alter- 
native constraints which can be put on the optimization 
are (i) an upper limit on the photon number, or (ii) a 
fixed mean photon number. The case where an upper 
limit N is put on the photon number was considered in 
Refs. [0, H , where the minimum phase variance was found 
to be 

V(4>) S3 2cN- p + (-z 1 )(2cp) 2/3 7V- 2(1+p)/3 , (4.1) 

where z\ is the first zero of the Airy function, and is equal 
to approximately —2.338. 

Here we will be considering the case of fixed mean pho- 
ton number. Let us take the phase to be zero and use 
the operator for the phase variance 

5 = 2- 2cos^ (4.2a) 

OO 

= 2- [l-ft(m)](|m)(m + l|+H.c). (4.2b) 

m=0 

In order to optimize for minimum phase variance with 
a fixed mean photon number, we use the method of un- 
determined multipliers. We have two constraints on the 
state: that the mean photon number is fixed and that 
the state is normalized. This then gives us an equation 
with two undetermined multipliers 



(5 + /iATM=HV>>. 



(4.3) 



We solve this as an eigenvalue equation for v with a fixed 
value of n, and the eigenstate corresponding to the mini- 
mum eigenvalue is the optimized state. The mean photon 
number can then be found from the state. The mean pho- 
ton number can be varied by varying /x, but cannot be 
easily predicted from u. 

This method was used by Summy and Pegg JTT| to 
find states with fixed mean photon number optimized for 
minimum intrinsic phase variance. They found that the 
minimum phase uncertainty was approximately 



((^0) } intrinsic 



c 



(n + e) 2 



(4.4) 



where C ps 1.88 and e ps 0.86. Therefore, the minimum 
intrinsic phase variance scales as n~ 2 . 

If the state is expressed in the number states basis 

00 

|V>=y>„|n), (4.5) 



then in terms of b n , Eq. fl4.3| ) becomes 

2b n - [1 - h(n)] (b n+1 + = {v - /in) b n . (4.6) 
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In the case of very large mean photon number, we can 
use a continuous approximation 



2h(x)y 



l dx 2 



[l-h(x)]^ = (v-nx)y, (4.7) 



where we are taking x = n and y{x) — b n . For large n, 
[1 - h(n)] ps 1, so 



fy 
dx 2 



+ y [2h{x) -i/ + ^]«0. 



Now define 

f(x) = 2h(x) — v + [ix. 
Then, expanding f(x) in a Taylor series around 
*o = (»/2cp)-^ +1 \ 



we find 

/(so) 
/'(so) 
/"(*<,) 
/'"(so) 



2h(x ) - v + fix , 
2h'(x ) + M = -2cpx/ _1 + = 0, 
2^"( a; o) = 2 C p(p+l)a ;o p - 2 , 
-2cp(p+l)(p + 2)x - p - 3 . 



(4.8) 
(4.9) 
(4.10) 

(4.11a) 
(4.11b) 
(4.11c) 
(4. lid) 



This technique requires that the number distribution has 
its maximum near xq, which will be justified below. 

Using the Taylor series for f(x), and defining f — 
f{x ), h = /^4o)/2, and f 3 = f'"(x )/6, the differen- 
tial equation ( |4.8| ) becomes 

-0 + [h{x - .T ) 2 + f 3 (x - x f]y w -f y. (4.12) 

Note that — fo = v — [2h(xo) + fJ-x p], s o the above equa- 
tion is equivalent to solving Eq. (4.3) as an eigenvalue 
equation for v with a fixed value of \i. Now Eq. ( 4.12j ) 
is equivalent to the time-independent Schrodinger equa- 
tion with energy eigenvalue E — — fp for a perturbed 
harmonic Hamiltonian H = Hp + Hi , where 



Ho — \fh 

Hi = be, 



d 2 



where 



Mp+i)] 1/4 ^" /2 , 



fl /A {x - Xq). 



The unperturbed solution is 



(4.13a) 
(4.13b) 

(4.14a) 
(4.14b) 

(4.15) 



where Hj are Hermite polynomials. This solution is only 
valid for \fj2x\ 3> 1, which requires that p < 2 in addi- 
tion to xq 3> 1. The energy eigenvalues are 



(4.16) 



The lowest energy eigenvalue Eq and corresponding 
eigenstate yo will minimize v and therefore minimize (S). 
From perturbation theory they can be expressed as 



2^ 



E 

fc=i 



yo 
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(fc|^ 3 |0) (o) 

— r— M . 



(4.17a) 
(4.17b) 



where \j) is the state corresponding to y^iO- Now it 
is easily shown using the properties of Hermite polyno- 
mials that the only nonzero terms are (1|£ 3 |0) = 3/2y2 
and (3|£ 3 |0) = V3/2. This then gives the lowest energy 
eigenvalue and eigenstate as 



E PS E, 



(0) 
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(4.18a) 
(4.18b) 



Now we can use these expressions to find the mean 
photon number as 



n=(x)= / yo(0(*o + / 2 1/4 02/o(£R 



X 



J>/2 



Vcp(p+ 1) 
As we are assuming p < 2, 



x 



1 



4i/cp(p + 1) 



,P/2-l 



(4.19a) 
(4.19b) 

(4.20) 



so that the mean photon number is close to xp 7 justifying 
the previous expansion around xp. Now we can find the 
minimum phase variance using 



£tn) r 



(4.21) 



/mm; 

and v = —fo + 2h{xp) + ^xp, /i = 2cpx Q p , and fo 
—Eq- We obtain the first two terms as 



2cn" p + A /cp(p+l)rr p/2 - 1 . 



(4.22) 



Note that the first term here is the same as the result 
when an upper limit is put on the photon number, but 
the second term scales as a different power of n. 

A particular case of interest is heterodyne detection, 
for which we find 



Whet ~ ^ 



1 



2n 3 / 2 ' 



(4.23) 



This is interesting because it differs radically from the 
result claimed by D'Ariano and Paris |l2j of 



2 \ min 
het 



1.00 ±0.02 
~1.30±0.02 ' 



(4.24) 



As the quoted errors suggest, this result was obtained 
entirely numerically in contrast to our analytical result. 
In Sec. VI we present our own numerical results and show 



that our analytical result is a far better fit than the power 
law of D'Ariano and Paris. 
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V. OPTIMIZED SQUEEZED STATES 

As an alternative to considering general optimized 
states, we can consider optimized squeezed states. There 
are three reasons for this: 

(1) Squeezed states are relatively easily generated in 
the laboratory, whereas there is no known way of pro- 
ducing general optimized states experimentally. 

(2) Squeezed states can be treated numerically far 
more easily than general optimized states. 



(3) It has been found numerically (see Sec. VI) that the 
phase uncertainties of optimized squeezed states are very 
close to those of optimized general states, and a partial 
theoretical explanation can be obtained by the following 
analysis. 

Squeezed states optimized for minimum intrinsic phase 
variance were previously considered by Collett [ [i3| , who 
found that the minimum phase uncertainty was approx- 
imately given by 



intrinsic 



ln(n) + A 

4^ ' 



(5.1) 



where A = | + 2 In 2 — \ ln(27r). The scaling of optimized 
squeezed states is therefore worse than the scaling of op- 
timized general states when we consider intrinsic phase 
variance. The difference is only a factor of Inn, however. 
We consider squeezed states of the form 



\a, C) — exp(aa^ 



a*a)exp[(C*a 2 



Co T2 )/2]|0). (5.2) 



Now we will take the phase to be 0, so a is real. Phase 
squeezed states have £ real and negative, so we can take £ 
to be re al. Then using the definition of phase uncertainty 
in Eq. (O) gives 



<(A^) 2 )=2-2^( a ,C|n)(n + l|a,C) 

CO 

+ 2j2h(n)(a,C\n)(n+l\a,0- (5.3) 

71=0 

In estimating the intrinsic phase uncertainty, Collett 
found the first two terms on the right-hand side to be 
approximately (n + l)/4ri 2 + 2 crfc(v / 2no), where no — 
ne 2 ^ . We therefore only have to determine an expression 
for the third term. Using h(n) w cn~ p as usual, the result 
we require is derived in the Appendix: 



Y^n- p {a,C\n){n + l\a,() « n~ p 



p(p+l) 
2n 



The phase uncertainty is therefore given by 



(5.4) 



((A0) 5 



no + 1 

An 



+ 2 erfc( % /2^) + 2cnr v 



pip + 1) 
2n 



Taking the derivative with respect to no gives 

d((Acf>) 2 ) _ 1 ^„-2n 



dn 



1 

i^ 2 



+ 2cn~ p 



irn Q 



pjp + 1) 
2n 2 Q 



(5.6) 

As the second term falls exponentially with no, it can be 
omitted. Then we find that the minimum phase variance 
occurs for 



n « 2^cp{p+l)n 1 ~P/ 2 



(5.7) 



Thus we find that no oc n 1 p / 2 , as was u sed in the Ap- 
pendix. Substituting this result into Eq. (5.5) gives 



((A0) 2 ) min » 2cn-P + ^/cp{p+l)fi- p / 2 - 1 . (5.8) 

We therefore obtain exactly the same first two terms for 
the phase uncertainty when considering squeezed states 
as we do when considering general states. 



VI. NUMERICAL RESULTS 



The analytic results from Sees. IV and have been 
verified numerically by calculating the optimized states 
for heterodyne measurements and adaptive mark I and II 
measurements. For moderate mean photon numbers the 
calculations were exact, except in that a cutoff at large 
photon numbers was used. 

For larger photon numbers an additional approxima- 
tion was that asymptotic expressions for h(m) were used. 
For heterodyne measurements there is an exact expres- 
sion for h(m) pL 



K c t{m) 



T (m + 3/2) 

(to + i) r (to + 2) ' 



(6.1) 



This form of the equation is very inconvenient for use in 
numerical work due to roundoff error. It is more conve- 
nient to use the asymptotic expansion. From the asymp- 
totic expansion of lnr(m) [H it is simple to find an 
expansion of /ihet(w) in powers of 1/ (to + 1) to as many 
terms as required. The first few terms are given by 



hhetim) 



1 



1 



8(m+l) 2 7 (to + 1) 5 



O (to 



(6.2) 



An expansion up to 12th order was used to determine 
hhetim). This expansion was found to be more accurate 
than using the formula directly for values of to>12. 

For adaptive mark I and adaptive mark II measure- 
ments, h(m) can be determined using methods discussed 
in Ref . |jj . For the mark I case hi (to) was determined ex- 
actly up to to = 3000. Further values were extrapolated 
by fitting an asymptotic expansion to the results below 
3000. The first three terms of the asymptotic expansion 
used were 



(5.5) hi(m) = 



l 



0.101562 0.0508 



8mV2 



,3/2 



+ 0(m~ 2 ). (6.3) 
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The uncertainties in these numbers were estimated us- 
ing various methods of fitting, and are indicated by the 
number of significant figures quoted. 

For the mark II case hn(m) was determined exactly 
up to m = 1000. Fitting techniques did not consis- 
tently give any higher-order terms than that obtained 
by the semiclassical theory, so for m > 1000 the formula 



hu(m) 



16 



m 



~ 3 / 2 was used. 



For very large mean photon numbers ^10 5 , it was not 
feasible to solve the exact eigenvalue problem, but an ap- 
proximate solution was obtained by using the continuous 
approximation of the eigenvalue problem and discrctizing 
it. In order to reduce the number of intervals required 
in the discretized equation, the equation was solved for 
three different numbers of intervals. The result for the 
continuous case was then estimated by projecting to zero 
step size, assuming the error is quadratic in the step size. 
For optimized squeezed states the full calculation was 
performed up to a mean photon number of about 10 6 , 
beyond which roundoff error became too severe. 

The results of the numerical calculation for the gen- 
eral optimized states are shown in Fig . [I], along with the 
analytical results obtained in Sec. IV and the power law 
of Eq. Q4.24D published by D'Ariano and Paris O. It is 



seen that the results for heterodyne measurements agree 
reasonably well with the power law of D'Ariano and Paris 
for moderate photon numbers (up to about 100). Above 
this, however, the agreement is extremely poor. This 
is presumably due to the fact that the numerical data 
used by D'Ariano and Paris seems to have been limited 
to maximum photon numbers only of order 100. In con- 
trast, our analytical result agrees very well for mean pho- 
ton numbers above about 100. The analytical result also 
agrees well with the numerical results for the adaptive 
mark I and II measurements. 

On the log- log plot of the phase variance it is extremely 
difficult to discriminate between different phase variances 
unless the difference is greater than a factor of 2. There- 
fore, in Fig. H we plot the parameter z, defined by 



z = «(A</>) 2 



2cn~ p )n p / 2+1 . 



(6.4) 



From the above analysis this parameter should converge 
to \J cp(p + 1) for large n. The z parameter is plotted for 
optimized general states and optimized squeezed states 
for the cases of mark II, heterodyne and mark I measure- 
ments in Figs. §(a), |](b), and|2|(c), respectively. 

In the adaptive mark II and heterodyne results in Fig. 
^(a) and ||(b), we can see that the results using the full 
eigenvalue solution match up very well with the contin- 
uous approximation results, demonstrating the accuracy 
of this technique. In addition, the squeezed-state results 
are extremely close to the general optimized state results, 
far closer than indicated by the first two terms derived 
above. Also note that in Fig. |^(a) the results do not 
agree closely (within 1%) with the asymptotic value un- 
til n > 10 8 , whereas the heterodyne results converge at 
a much lower photon number, with good agreement for 
n«10 4 . 



10" 



10 



E10~ 




FIG. 1: Plot of the phase variance for general optimized states 
via complete eigenvalue solution (crosses) and continuous ap- 
proximation (circles) vs mean photon number n. The asymp- 
totic analytical expression from Sec. [Ey| is also plotted (con- 
tinuous lines). The three phase detection schemes are (a) 
mark II adaptive, (b) heterodyne, and (c) mark I adaptive. 
The power law claimed by D'Ariano and Paris for heterodyne 
detection is also plotted in (b) (dash-dot line). 
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FIG. 2: Plot of z parameter from Eq. ( p.4[ ) for gen- 
eral optimized states via the complete eigenvalue solution 
(crosses) and continuous approximation (circles), and opti- 
mized squeezed states (dashed lines) , vs mean photon number 
n. The theoretical asymptotic value of \J cp(p + 1) is shown 
as the continuous horizontal lines. The results shown are for 
mark II adaptive phase measurements in (a) , heterodyne mea- 
surements in (b), and mark I adaptive measurements in (c). 
The predicted z for the mark I case taking into account the 
second term for hi(m) is also plotted in (c) (dotted line). 



In the adaptive mark I results in Fig. |](c), there is 
again good agreement between the squeezed-state results 
and the general optimized state results. However these 
results do not approach the asymptotic value at all in 
this case. The reason for this is that there is a higher- 
order term in hi(m) which is of order m , as shown in 
Eq. (S.3). This term is of lower order than the second 
term in the expansion Q4.22D of the phase variance, which 
is of order to -5 / 4 . If the higher-order term in h\(m) is 
taken into account in predicting the asymptotic value, 
then there is good agreement with the numerical results. 

Although the adaptive mark I scheme has the poor- 
est performance for large mean photon numbers, note 
from Fig. |l| that for small mean photon numbers (of order 
unity) it is actually the best scheme, having the smallest 
minimum phase variance. This is to be expected from the 
results of previous work |l], showing that for a maxi- 
mum photon number of 1, the adaptive mark I scheme is 
actually the best possible. 

The significance of these results is, first, that the two 
terms given by the theory for optimized general and 
squeezed states are correct, and, second, that the phase 
uncertainties of optimized squeezed states are extremely 
close to those for optimized general states. This result is 
of great importance, as it means that in numerical and 
experimental work squeezed states can be used rather 
than general states. 

This is an advantage in numerical work because only 
the squeezing parameters need be considered, rather than 
the entire state. This means that, for example, a numer- 
ical evaluation of different feedback schemes is feasible, 
and this work is being carried out now. In experimental 
work it is an advantage because squeezed states can be 
produced experimentally, whereas arbitrary states can- 
not. This means that it is possible to produce states 
experimentally that are very close to optimized for the 
different measurement schemes. 



VII. CONCLUSION 

We have derived asymptotic analytical expressions for 
the minimum phase variances obtainable under vari- 
ous detection schemes, for states constrained by their 
mean photon number. The detection schemes consid- 
ered were heterodyne detection (a standard scheme) , and 
two single-shot adaptive schemes (first proposed in Rcfs. 

[|), called mark I and mark II. Numerical results 
confirm the correctness of the first two terms in the 
asymptotic expansion, except in one case (adaptive mark 
I), where the second term was not expected to be cor- 
rect. Furthermore, analytical and numerical results show 
that essentially the same results may be obtained using 
squeezed states, rather than completely general states. 
This is an important result from both theoretical and 
experimental standpoints. 

As expected, the minimum phase variance for adap- 
tive mark II measurements was much smaller than that 
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for the standard technique of heterodyne detection, for 
large mean photon number n. In particular, the leading 
term in the former scaled as n -3 / 2 compared to YT X in 
the latter. The claim by D'Ariano and Paris |0] that the 
heterodyne phase variance scaled as n -130 (or, as stated 
in their abstract, n~ 4 / 3 ) was proven wrong. This rein- 
forces the position of adaptive mark II phase measure- 
ments as the best known phase measurement scheme. 

In Ref. JjJ it was shown that a lower bound for 
the phase uncertainty introduced by measurements is 
(lnn)/(4n 2 ). For large photon numbers this is a lot 
smaller than the phase variance of l/(8n 3 / 2 ) introduced 
by mark II measurements. It therefore may be possible 
to obtain a higher power in the scaling law by modifying 
the measurement method. The most promising modifi- 
cation is using a different feedback phase, and this is cur- 
rently under investigation numerically via the solution of 
stochastic Schrodinger equations |l(], [l5) . This is possible 
even with very large photon numbers if one uses squeezed 
states, because these remain squeezed states even under 
the stochastic evolution. The results obtained in this pa- 
per justify this approach, as the variances obtained for 
the cases of general states and squeezed states were al- 
most indistinguishable. 
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Rearranging this and squaring gives 



APPENDIX: DERIVING EQ. ( \5.4\ ) 

We wish to evaluate the following sum: 

oo 

J2Hn)(a,C\n)(n + l\a,0- (Al) 
n=0 

We can do this in the following way. The number state 
representation of squeezed states is given by jl6| 

(n\a,() = (nlv)- 1/2 (v/^) n/2 H n {m^r 1/2 } 

xexp[-|/3|72 + (^72M)/? 2 ], (A2) 

where 

[i = coshr, v = e 1 ^ sinhr, (3 — u[i + a* v. (A3) 

Here r and </> are the magnitude and phase, respectively, 
of while H n are Hermite polynomials and satisfy the 
recursion relation Hi 



H n+ i{x) - 2xH n (x) + 2nH n - l {x) = 0. 



(A4) 



This means that the number representation of squeezed 
states satisfies the recursion relation 



(n + l\a, C)/iVn + 1 - (n\a, 0/3 + (n - l\a. ()v-s/n = 0. 

(A5) 



2 o2 



|(n + l|«,C>|V>+l) = IHa,C)r/3 

+ \(n - l\a, C}| Vn - 2(a, fln) (n - l|a, ()/3v^ 



(A6) 



This expression is only true for real squeezing parame- 
ters. Multiplying this by n k and summing gives 



2(3vJ2(n + l) k + 1 / 2 (a,(\n){n + l\a,0 

n=0 

= I3 2 {n k ) + v 2 {{n + l) k+1 ) - fi 2 (n(n - l) k ). (A7) 

Now let us take — p = k + i. In this case we find that 
the shift of indices cannot be performed exactly, but the 
contribution from terms near n = will be negligible. 
Also some of the terms above diverge near n = 0; how- 
ever, the divergent terms are the extra terms produced 
by the shift of indices, and in the following expansions 
the behavior near n = is ignored. 

Taking — p = k + ^ and considering the deviation from 
the mean photon number gives 



2ppJ2(n+l)-i>{(3\n){n+l\ 



n=0 



w (3 2 ((n + Any^+W) +is 2 ((n+l + Any^ 1 ^) 
-^((n-l + An)-^- 1 ^) 

-Li 2 {(n- 1 + An)- (p+ W). (A8) 
Expanding this in a series in An gives 

oo 

2/3i^(n+l)- p (/3|n)(n+l|/?) 



n=0 



3=0 



tm^j r (p+j-i/2)! ( ^ 2 __ fa+ , +1/21 



(p- 1/2)1 



Ml > )+ (p- 3/2)! 

x \v 2 {n + l)-(f+J-l/2) _ M 2 (n _ ^-(p+i-l/2)- 



(A9) 



N ow we have an expression we can use to evaluate Eq. 
([Al]). Recall that for generalized measurements we have 
the asymptotic expression /i(n)«cn~ p . This is equivalent 
to h(n) ~ c(n+ l) _p , as the difference is of higher order. 
It is easily shown that for squeezed states (An 2 ) = a 2 (fj,— 
v) 2 + 2ji 2 v 2 . Therefore, using the first three terms of Eq. 
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( |A9j ) above, we find 

oo 

2/3uJ2in+ l)- p (a,(\n)(n + l\a,Q 

71=0 

«/3 2 n"(P+ 1 /2)_ (U 2 ( -_ 1) -( P +i/2) 
+ v 2 {n + 1)-(p-V2) _ tffr _ i)-(p-i/2) 
"a 2 ^-^) 2 



+ 



, 2 2 
+ fJ, V 



P 



2P+1 



x (p 2 n-^+ 5 ^ - /i 2 (n- l)~(P+5/2)) + (V _ 1^ 
x[^ 2 (n+ 1)-(p+ 3 / 2 ) - M 2 (n - 1)-(p+ 3 / 2 )]|. 

(A10) 

At this stage the main problem is to determine which 
terms should be kept. This depends on how no scales 
with n. If the state is optimized for minimum intrinsic 
phase uncertainty, then uq oc ln(n) fu| . If we carry the 
derivation through using this result to estimate the order 
of the terms, then we obtain the result no oc n 1-p / 2 . 
If we use this to estimate the order of the terms, and 
omit all terms on the rig ht-ha nd side of order higher than 
--(p-i/2)^ then Eq. ( |aTo|) simplifies to 



DO 



n=0 

1 



2n 



p(p+l) 



(All) 



The first term on the right-hand side is not of higher or- 
der than fi - ^" 1 / 2 ) /no if p < 1. If we were estimating the 
order from the parameters optimized for minimum intrin- 
sic phase variance then the first term would be omitted. 
Now we can expand —2/3v to give 



-2j3v 



If we were estimating the order from the parameters op- 
timized for minimum intrinsic phase variance, the term 
no/n would be omitted. This then gives us 



l-^)£(n + l)->,C|n)<n + lhC) 
n / i — ' 



-non 



orio 2no 



P(P+1) - J 



(A13) 



The two terms that would be omitted if we were consid- 
ering parameters optimized for minimum intrinsic phase 
variance just cancel, giving the simple result 



J2( n + i y P (uX\n)(n + l\a,()~n-' p 



71=0 



2no 
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